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Abstract 

" 1 

^ ■ We determine exact recurrence relations which help in the evalu- 

ation of matrix elements of powers of the radial coordinate between 
Dirac relativistic hydrogenic eigenstates. The power A can be any com- 
plex number as long as the corresponding term vanishes faster than 
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r as r — > oo. These formulas allow determining recursively any ma- 
trix element of radial powers — r A or /3r A , j3 is a Dirac matrix — in 
terms of the two previous consecutive elements. The results are useful 
in relativistic atomic calculations. 

1 Introduction 

Matrix elements of the radial coordinate have been very important since the 
early days of quantum mechanics [TJ|21Ellll0inilZllH]- These quantities 
are a crucial link between theoretical predictions and the observed facts. For 
example, they are needed for studies of the behaviour of Rydberg electrons 
in external fields, or in any calculation involving the multipolar expansion of 
electromagnetic fields 0. On the other hand, the interest in calculating rel- 
ativistic effects on high Z atomic or ionic transitions [10 i requires calculating 
matrix elements between relativistic eigenstates of the Dirac hydrogen atom. 
Experiments with precisions of the order of 10~ 3 eV or less are now standard. 
For example, experiments performed using new photon sources for the Opac- 
ity and Iron Projects, or experiments with merged beam techniques, or the 
use of specialized laser sources jTTJ H21 EH] . All of these mean that relativistic 
effects can be quite easily observed in atomic, molecular or ionic processes 
and, therefore, that techniques for evaluating relativistic expectation values 
are quite useful — as has been known for some time now [14j. 

The nonrelativistic evaluation of matrix elements has a long and succesful 
history, but the relativistic efforts have mainly come from the last 35 years 
or so. We have to pinpoint however, that even these relativistic calcula- 
tions, even exact ones, are an approximation, for the only way of calculating 
exactly relativistic effects in the interaction of atoms with electromagnetic 
fields is using the full QED formalism. See, for example, ^S]- But even so, 
many calculations use a nonrelativistic approach recurring to the Schrodinger 
equation — which is also a valid approximation ^2]- Even in the relativistic 
quantum mechanics approximation, in many circunstances it is much better 
to have recursion relations between power-of-r terms than to have to deal 
with the cumbersome and complex formulas that stem from the exact eval- 
uation of the matrix elements ^7j. Furthermore, all what is needed in many 
ocassions is to have the matrix elements for exponents in a certain range, so 
it is very convenient to have recursion formulas ^H]- For instance, the be- 
haviour of Rydberg electrons in external fields, or the long range interactions 
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of ionic cores ^H]- Thus, we have been investigating such formulas in Dirac 
relativistic quantum mechanics. We have discovered various of such expres- 
sions since, perhaps not surprisingly, there appears to be more independent 
recursions. We have been applying an hypervirial inspired technique that 
have produced excelent results in nonrelativistic quantum mechanics [211] to 
uncover some of the recurrence relations, but we are now combining that 
technique with operator algebra [21j to generalize the previously obtained 
recurrence relations to more compact recursions; equations (jH)) and (|45|) in 
section [31 An advantage of these recursions is that they are able to relate 
only three consecutive powers of r or (3r — to be compared with the results 
in 0HH] 

Our paper is organized as follows. In section [21 the relativistic quantum 
mechanics of an hydrogen atom is reviewed since all the discussions that 
follows use its energy eigenfunctions and eigenvalues. A useful feature of 
our discussion is that we express compactly the relativistic eigenfunctions in 
terms of generalized Laguerre polynomials of noninteger index [221 EH] ■ The 
relativistic recurrence relations of matrix elements of powers of the radial 
coordinate are derived in section [21 Section 0] gathers the definitions of the 
symbols we use for writing compactly these recurrence relations. Section [3] 
contains our conclusions. 



2 The relativistic hydrogen atom 

The radial wave function of an electron in an hydrogen atom is 



iG nje (r) 



(1) 



where F n j e (r), G n j e (r) are respectively called the big and the small compo- 
nents of the spinorial wavefunctions n — 1, 2, 3, . . ., j = 1/2, 3/2, 5/2, . . . 
is the total (orbital plus spin) angular momentum quantum number of the 
electron, and e = (— 1) J+ ^ 1//2 , I — j ± 1 is the orbital angular momentum 
quantum number, the sign is chosen according to whether I refers to the big 
(+) or to the small (— ) component. This wave function is a solution of the 
Coulomb radial Dirac equation 



ca r 



hf3 ( 1 
p r -i — e I j + - 
r V 2 



\[>(r) = E^(r) 



(2) 
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where Z is the atomic number, m the electron mass, q e is the magnitude of 
the electron charge, and the subscript r refers to the radial part of the wave 
function. The Dirac matrices are 



a. 



-1 
-1 



1 
-1 



and the radial momentum operator is 



p r = -in - + — . 
\r dr 



(3) 



(4) 



The components of the eigenfunctions, F and G, are the solutions of the 
following differential equations [221 EH 125] 



and 



v dp P J 



-v + 



Zap 



F{P) 



(5) 



where 



\ dp p 



1 



Hp) 



1 Zap 

y p , 



G{P) 



p = kr, k = —\/m 2 c A — E 2 , 
nc 



E 



v 



(6) 



(7) 



mc 2 + E 

The solutions of these coupled equations — that can be expressed in terms 
of the Sonine polynomials, in terms of Laguerre polynomials, of non-integer 
index, or as special cases of the hypergeometric function, see |221 EH I2H E3 
EE! f° r details — can be written as 



F(p) = +CVmc 2 + Ep s exp(-p) A n L 2 n s (2p) + B^P) 



and 



G(p) = -CVmc 2 -Ep s exp(-p) A n L 2 n s (2p) - B 2 n s _^P) 



(8) 



(9) 



where the L"^(x) are Laguerre polynomials of noninteger index, the numbers 
A n and B n are 



+ 2) +s + 



Zcti 



n 



Za.} 



V 



1/2 



(10) 
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B r , 



(n + 2s) 



ZctF 

s h n 



Zai 



v 



1/2 



where s 



(j + 1/2) 2 — Z 2 a 2 F . The normalization constant C is 

h2 s - 1 I 



C 



n\k 



Za F c 2 V 2m 3 r(n + 2s + 1) 
The energy levels are given by the usual expression j2H |27] 



(11) 



(12) 



E 



E, 



1 + 



Z 2 a\ 



n 



j-l/2 + J(j + l/2) 2 -Z 2 a F , 



-1/2 



(13) 



where n = 1,2,3,..., op = q 2 /hc ~ 1/137 is the fine structure constant. 
It is convenient to point out that the quantum number e we use is related 
to the often used k ^Hj as k = —e(j + 1/2). The quantum number k is an 
eigenvalue of the following operator 



K 



(3 (h 2 + s 



•L 



where X = 2S = h 



(14) 



and a = (a x ,a y ,a z ) is the standard 3-vector of the Pauli matrices. For the 
sake of simplicity we often use the single symbol a to mean all the three 
quantum numbers n a ,j a ,e a , and write the radial wave functions in ket (or 
bra) form as ^f a = \a). 

We calculate recursions beween matrix elements of the form 



(a'\(3 b r x \a) = J dr r 2 ^l{r)(3 b r x ^ a (r) 



(15) 



where W a , = (F*, /r, —iG*,/r), b = 0, or 1, and A is a possibly complex 
exponent. The recurrence relations (|2l?j) . (jSSJ), (JSHJ) -(JISJ) hold true if S1 + S2 + 
1 + |A| > dZ|EE| where as before the s b = +^J(j b + 1/2) 2 - Z 2 a 2 F , 6 = 1,2 
are real numbers; this conditions is basically that any integrand goes to zero 
faster than 1/r as r — > 00 f£F\ . 

One of the advantages of the approach developed here is that is not only 
applicable to Dirac hydrogenic wavefunctions but can also be used within 
the quantum defect approximation, as discussed by Owono Owono, Kwato 
Njock, and Oumarou jH], see also |7|. This feature is important since the 
quantum deffect approximation describes accurately the behaviour of Ryd- 
berg electrons in atomic systems [3 IHl ESI EDI EI] • 
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3 Relativistic recursion relations 



Using hypervirial methods we have recently obtained a useful form of re- 
currence relations for atomic eigenstates. In this work we want to discuss 
the additional use of operator algebra to obtain alternative and indepen- 
dent expressions for recurrence relations of atomic expectation values. Let 
us begin with the radial Dirac equation for an electron in an arbitarry radial 
potential V(r), where the index a serves to distinguish Hamiltonians with 
different quantum numbers and it is convenient to distinguish wavefunctions 
with different energy eigenvalues. 

Let us begin with the Dirac radial equation written with a certain radial 
potential that afterwards will be identified with Coulomb's 



ca r 



hp ( 1 



+ (3m a c 2 + V a (r)j * a (r) = E a ^ a (r), 

(16) 

here H a is the Dirac Hamiltonian, the energy eigenvalue E a is given by (|13|). 
and a is an index that will be quite useful in what follows. 

To begin with, let us consider an arbitrary radial function f(r) and define 
£(r) = H 2 f — fH\. To obtain the recurrence relations, let us establish that 



2r \ ar r J 2r 

'A 



+ ^ 2 (^r) f+(c^m- + v-)f 



ifica r 



\ ar r 2r 2r 



(17) 



where, A a = e a (2j a + 1) and, if X is any symbol of interest, we are — and 
will be — using 



X ± = X 2 ±X 1 . (18) 

The following five identities are not too difficult to establish by the use 
of a little operator algebra 
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d 



J * 



H 2 (-ica r f) + {-icorf)H x = -c 2 h [2f— + f + 2 J - - — (5f 



dr 



r 2r 



+ ica r (c 2 (3m — V + ^j f, 



(19) 



H 2 f - fH x = -ihca r If + —f3f) + (c 2 (3m- + V~) f, (20) 



d 



2r 

J , A" 



Htf + fHi = -i%ca r hf- + f + 2 J - + —Pf\ + (c 2 (3m + + V + ) f, (21) 



H 2 V-f-V~fHi = -ihca r [V-f + (V-)'f + —V-(3fj+(c 2 p m - + V~) V~f, 

(22) 



and 



H 2 (— ica r f3— j + (— ica r (3— j ifi 



-e 2 ft 



\ r r 2 / 2r r 



J 



(23) 



To get the relation we are after, we need to eliminate the terms involving 
ica r and the terms involving d/dr of equation (|17jl. To this end, from (fTT)|) 
we extract the term in the left hand side (LHS) of the following equation 



.A- 



f 



- c 2 h 2 p—\2f- + f + 2 J - 
2r \ dr r 



'A 



hca r f3— (c 2 pm- - V + ) f 



A 



c h hH / - H 2 {-ihca r (3—f) + {-i%ca r P—f)H x . (24) 



2r 



2r 



2r 



Furthermore, from (|20p . we extract 



- z/ica r ( /' + ) = - (#2/ - fH x ) - (c 2 Pm~ + /, (25) 



2r 



from (|21j) . we extract 
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d 



iha r c 6 (5m- 2f— + f + 2- 



dr 



H 2 c 2 m-{3f + c 2 m-{3fH 1 



A 



+ %hca r c 2 mT / - f c 2 (3m + + V + ) (3f } 

2r '• 



from 



we extract 



(26) 



and from 



- ihca r [y-f + (V-)'f) = [H 2 V-f - V-fHi 

+ ihca r —V-(3f - (c 2 f3m- + V~) V~f, 
2r v / 



we extract 



(27) 



( A- f 

H 2 —ic%a r f3— 

\ 2r r 



A" A 

—icha r [3— H\ 

2r r I 



icha r — (c 2 (3m~ - V + ) (3f = -c 2 h 2 — 
2 r \ I 2r 



( r) 2r ^ 



.(28) 



The LHS of equations (|2^jl to (jSHj) are to be substituted into equation (fTTj) . 
Then, after introducing explicitly the Coulomb potential V\ = V 2 = —Zq^/r, 
imposing the equality of the masses mi = m 2 (since H\ and H 2 describe 
the same system), substituting the arbitrary function f(r) for the potential 
function r A , and taking matrix elements between states (2| and |1) at the 
end, these steps yield the recurrence relation 



(E 2 - E,) 2 (2\r x \l) + k 2 (2\r*- 2 \l) = l (2\(3r'\l) + l 2 (2\(3r*- 2 \l) , (29) 
where 



h 



c 2 h 2 



c 2 h 2 X (A - 1) 

-2c 2 m(E 2 - E x ), 
c 2 h 2 



4 



2A~ + A A + - A 



(30) 
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Now, following a strictly similar procedure but substituting equations 
(Hi, PI, (J23J), and the following one 



H 2 V~f + V-fHi = (c 2 pm + + V + ) V~f 
— ihca 



2V~f£- + V~f + (V-)'f + 2V- f - + V-^-Pf 
ar r 2r 



(31) 



into the next equation 



\ ar r 2r r , 



.2*2 



A+ A" 



2r 2r 



/ + (c 2 /3m+ + V + ) (c 2 Pm- + V~) f 



ihca r 



f + ^-Pf) (y + - c 2 Pm-) + c 2 (3m-f + c 2 m -^-f 



2V-fj- + V~f + (V-)'f + 2V~ f - + V-^Pf 
ar r 2r 



(32) 



again introducing explicitly the Coulomb potential V% — V% — —Zq 2 /r, mak- 
ing the masses equal mi = m 2 , substituting the potential function r A , and, 
at the end, taking the matrix elements, we obtain the following recurrence 
relation 



(E 2 + E 1 ) (E 2 - E 1 ) (2|r A |l)+M 1 (2|r A - 1 |l)+M 2 (2|r A - 2 |l) = N 2 {2\(3r x - 2 \l), 

(33) 

where 



Mi 


= 2E~Zq 2 E , 




M 2 


= c 2 h 2 \{\-2) 


A+A" 

4 


N 2 


c 2 n 2 

2 


A" + A 


(A+ - A-)] 



(34) 



The new recursions [equations (}2*9*j) and ([33)1] can be more useful if the matrix 
elements of terms r A were uncoupled from the matrix elements of terms /3r A . 
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To disentangle such relations, we need the following three previously reported 
recursions [UJ [TH] , 



c (2|r A |l) = c 1 (2|r A - 1 |l)+c 2 (2|r A - 2 |l)+c 3 (2|r A - 3 |l) = d 2 (2\f3r x - 2 \ l)+d 3 (2\/3r 

(35) 

e (2|/3r A |l) = 6 (2|r A |l) + 6 2 (2|r A - 2 |l) + e 1 (2|/?r A - 1 |l)+e 2 (2|/5r A - 2 |l), (36) 
and 

^ 2 (2|r A - 2 |l) =p 2 (2|/3r A - 2 |l)+p 3 (2|/3r A - 3 |l). (37) 

where the coefficients b, c, d, e, g, and p are all defined in [TH] . 

Using now and (|37jl in ()29|). (|35|) and (J3BJ) . we get the following six 
uncoupled recurrence relations, 

A (2|r A |l) = A 1 (2|r A - 1 |l) + A 2 (2|r A - 2 |l> + A 3 (2|r A - 3 |l>, (38) 
B (2|r A |l) = B 1 (2|r A - 1 |l) + B 2 (2|r A - 2 |l) + B 3 <2 |r A ~ 3 | 1 ) . (39) 

C (2|r A |l) = C 1 (2|r A - 1 |l) + C 2 (2|r A - 2 |l> + C 3 (2|r A - 3 |l> + C 4 (2|r A - 4 |l), (40) 

D (2|/5r A |l) = D 1 (2|/3r A - 1 |l)+D 2 (2|/3r A - 2 |l) + D 3 (2|/3r A - 3 |l)+D 4 (2|/?r A - 4 |l), 

(41) 

E (2|/5r A |l) = E 1 (2|/5r A - 1 |l) + E 2 (2|/3r A - 2 |l> + E 3 (2|/3r A - 3 |l>, (42) 

and 

F (2|/5r A |l) = F 1 (2|/3r A - 1 |l) + F 2 (2|/3r A - 2 |l) + F 3 (2|/3r A - 3 |l>. (43) 

The equation obtained from substituting (J57)l into (JSfiJl [equation (jj2| above] 
coincides with one of the previously obtained recursions [equation (18) in Ref. 
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18j]. This is nice since it serves as a way of checking the above recursions. 
The coefficients in equations (J2ED, ©, (001), (ED), (H2J), and (USD are all 
gathered in section HJ 

The equations ()38|) and (J3TJJ) above express two forms of a recurrence 
relation relating matrix elements of r A with matrix elements of r x ~\ i = 
1, 2, 3. This is a manifestation of the richer behaviour - - though at the 
same time more restricted — of relativistic quantum mechanics as compared 
with the nonrelativistic theory. We can combine anyway equation (J3*%J) with 
equation (J3TIJ) to get the simpler relation 

A similar situation occurs with equations and above; they can be 
combined to yield the simpler relation 




(45) 

Equations (jjHj) and (JUJ) are simpler, and thus potentially more useful, than 
those reported in [01 HE] • 



4 The coefficients in the recursion relations 

The coeficients used in the recursion relations [equations (J38|) . (JHHJ), (jH3J), 
(JUJ), (J12J, and (JlHJl ] are defined as — the various symbols used to write these 
equations are defined at the end in equation (JS2J) - 



A 
Ai 

A 2 
A 3 



hE+(E~) 2 A- hE+E-A~W 



D 

2E~K 



R 



D 



(X-l)K 



2hE+A- 



h c 2 hw 



D + Ahmc 2 [R - h(A+ - A- 



4R ^ D + 4hmc 2 



+ 



AhZq 2 A- 



c 2 hA-p 



+ 



Q 



D + Ahmc 2 



R - h(A+ - A" 
g + 12^ 2 (l-A) 
2h 2 [A+ - A(A+ - A-)] 



(46) 
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The A a are the explicit coefficients in equation ()38j). 



B = -8Zq 2 e D E L 



B 1 = AhXF 



c 2 KR ' 

lQ(Zq 2 e ) 2 E-D A+E+E-L A(E + ) 2 E-D 



c 2 ft [A+ + A(A+ + A - )] R c 2 h[R + 2ft(A+ - A - )] ' 

B 2 = -2A£D-. 



, 4(A 2 - 1) - A+A- 2Zq 2 e E-A+L 2Zq 2 e E-E+D 



[A+ + A(A+- A-)] R c 2 %[R + 2ft(A+ - A~ ' 

B 3 = -c^(i_A)L-— A + L(S-Ah\)-2E + D 9r t~t^ \7> (47) 

v ; AR v ; 2ft [i? + 2ft(A+ - A - )] V ; 

The B a are the explicit coefficients in equation (J39*j) . 



4mE+ (£ 

C — 



\2 



Ci 



ft[# + 2ft(A+- A-)]' 

8mZg e 2 (£") 2 
ftp? + 2ft(A+ - A - )] ' 



.ii'/ 

C ? = -mc l E 



S + 12ft A 



ft[_R + 2ft(A+ - A" 

3 = 2Zq 2 e E- 1 



2 „_R + hA- 



R ' 

C 4 = i_(,S-4ft 2 A) (.R + ftA-) -c 2 ft 2 A(A-l) + ^-A+A-. (48) 
4ri \ / v / 4 

The C a are the explicit coefficients in equation (J3Uj). 



D \A+ J 

D \A+ D V J ' 



D 2 = -|(^) 2 (z 2 ) 2 A^ + ^)c 2 (T + ft)P+^A-^ 
D 3 = ^Zq 2 A-UP--J—(T-2h)Y+ Q 



2 ^ e hA+D v ' D + Ahmc 2 ' 

Zq 2 YE + (E-f 

° 4 ~ D + 4ftmc 2 ' ( 9) 
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The D a are the explicit coefficients in equation (|41j) . 



E = 2E+D-^F(T + h), 
At, 

E t = -AZq 2 e D + — Zq 2 e XFU, 

c 2 Ti 2c 2 
E 2 = _£^A+L-^L(l-A)(T-n) 

Es = -^(l-A)LC/, 
The E a are the explicit coefficients in equation 



4mc 2 /_ \3 . 

Fo = (r + R), 

F t = Zq 2 e (E-) 2 U, 

r 2 F + r 2 r 2 h 

F 3 = jZ 9e 2 SC 

The F a are the explicit coefficients in equation (|43|). 

The above equations are written in terms of the following symbols. 



D = hA~E~ - Amc 2 Ti\ 

T = A+ A + f hE + + h\, 

F = (E-) 2 -4m 2 c 4 , 

K = Zq^hE-A-, 

L = Ah 2 x 2 -h 2 (A-) 2 , 

U = ^(A+ + A-), 

W 7 = (1-A) AB+A+ 



A-E- -4mc 2 A' 
F = 2c 2 h 3 Zq 2 e (A - 1) (A - 2) A" 
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Y 



hA+ 



P 



D 2 



Q 

s 



R 



c 2 h 3 Zq 2 e (X- I) A~A + , 
Ah 2 \{\-1) -h 2 A + A-, 
HA- + Xh (A+ - A~) . 



(52) 



The coefficients appearing in equations and are respectively 
expressed in terms of definitions (|4l)]). (jUj) , and (|3U)l. (f5Tj) . 

5 Conclusion 

In this paper we have derived exact recurrence relations between general 
non- necessarily diagonal matrix elements of powers of the radial coordinate. 
These recursions relate any three consecutive powers of r or of j3r. The states 
used for evaluating the matrix elements are radial completely relativistic hy- 
drogenic eigenstates. The derivation was done employing a technique inspired 
in the hypervirial method and using some operator algebra [21 EH IU [2D] • The 
matrix elements analysed here are all gauge invariant [8 . The relations ob- 
tained may have different uses in many interesting atomic calculations, as in 
the calculation of transitions between Rydberg states beyond the semiclassi- 
cal approximation [HJ E2] . 
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